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ABSTRACT 

The linear dispersion relation for the magnetorotational instability (MRI) is derived including 
finite Larmor radius (FLR) effects. In particular, the Braginskii form of the ion gyroviscosity, 
which represents the first-order FLR corrections to the two-fluid equations, is retained. It is 
shown that FLR effects are the most important effects in the limit of weak magnetic fields, 
and are much more important than the Hall effect when Pi 3> 1, where /3, is the ratio of the 
ion thermal pressure to the magnetic pressure. FLR effects may completely stabilize even MRI 
modes having wavelengths much greater than the ion Larmor radius. Some implications for 
astrophysical accretion disks are discussed. 

Subject headings: accretion, accretion disks — instabilities — MHD — plasmas 



1. Introduction 

The magnetorotational instability is a local in- 
stability which may be present in accretion disks 
having sheared azimuthal flow and a weak mag- 
netic field (Balbus & Hawley 1991). Turbulence 
resulting from this instability is thought to play 
an important role in the radial transport of angu- 
lar momentum in such systems (Balbus & Hawley 
1998). Here we explore the modifications to the 
local, two-fluid theory of the MRI in the linear 
regime due to finite Larmor radius (FLR) effects. 

Rosenbluth et al. (1962) have shown using ki- 
netic theory that FLR effects can be stabilizing to 
"weakly unstable" modes — defined as modes hav- 
ing a linear growth rate much smaller than the ion 
cyclotron frequency — even when the mode wave- 
length is much larger than the ion Larmor radius. 
It was later shown by Roberts & Taylor (1962) 
that this result could be obtained from fluid theory 
by retaining the gyroviscous stress component of 
the ion pressure tensor. In most physical scenarios 
the MRI is weakly unstable in the sense of Rosen- 
bluth et al., and indeed we show that the MRI may 
be completely stabilized by gyroviscous effects at 
scales much larger than the ion Larmor radius. 
In some cases this stabilization significantly con- 
strains the spectrum of linearly unstable modes. 

The effect of the Hall term, which accounts for 



differences between the electron and ion fluid ve- 
locities, has been examined previously by Wardle 
(1999), Balbus & Terquem (2001), Salmeron & 
Wardle (2003), and Krolik & Zweibel (2006). In 
particular, it was found that the Hall effect may 
be either stabilizing or destabilizing, depending on 
whether the equilibrium magnetic field is aligned 
or anti-aligned to the equilibrium angular veloc- 
ity. It was also found that the Hall effect is im- 
portant only when the ion cyclotron frequency is 
comparable to, or smaller than, the orbital fre- 
quency. This situation may occur in early galaxy 
formation where the magnetic fields are still weak, 
or in weakly ionized protostellar disks. Krolik & 
Zweibel (2006) have suggested that, in this limit, 
short-wavelength modes are likely suppressed by 
viscous or resistive damping, leaving only slowly 
growing, long-wavelength modes as the magnetic 
field get sufficiently weak. However, their analy- 
sis is restricted to low-/^ plasmas as they do not 
consider FLR effects, which we show to be much 
more important than the Hall effect in the weak- 
field limit. The strong FLR stabilization of the 
MRI in the weak-field limit may have important 
implications for the possible role of the MRI in the 
amplification of weak, primordial magnetic fields. 

The gyroviscous stress is defined as the trace- 
less, perpendicular part of the ion stress tensor 
which does not depend explicitly on the collision 
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frequency (Ramos 2005). In typical cases, this 
stress arises primarily from variations in particle 
drift velocities over the scale of a Larmor orbit 
(Kaufman 1960). However, other effects may con- 
tribute to this stress, including gradients in heat 
fluxes. Braginskii's form of the ion gyroviscous 
stress is appropriate for collisional plasmas (in the 
sense that the ion mean-free-path is small com- 
pared to the hydrodynamic perturbation length- 
scale), and in the limit where the ion cyclotron 
period is short compared to collisional and hy- 
drodynamic time-scales (Braginskii 1965). More 
general, but more complicated, expressions for the 
gyroviscous force have been derived, which are ap- 
plicable to a broader range of collisionality regimes 
and dynamical time scales; Ramos (2005) and ref- 
erences therein provide derivations and discussions 
of these alternate forms. We choose to work with 
Braginskii's form here because it is (relatively) 
simple and applies to a broad range of astrophys- 
ical objects. 

The MRI in the collisionlcss regime, where 
the collisional mcan-free-path is greater than the 
mode wavelength, has been explored by Quataert 
et al. (2002) and Sharma et al. (2003) using kinetic 
closures. Islam & Balbus (2005) have extended 
the single-fluid MHD treatment to lower collision- 
ality regimes by including the Braginskii (1965) 
form of the parallel viscosity, and have obtained 
results similar to those obtained using kinetic clo- 
sure. These various analyses have found the linear 
growth rate to be enhanced by a factor of order 
unity at lower collisionality when an azimuthal 
component of the magnetic field is present, but 
the criterion for instability was found not to differ 
from the MHD result. We do include the paral- 
lel viscous stress for completeness in our analy- 
sis, as formally it may be larger than the gyro- 
viscous stress. However, for the sake of simplicity 
we restrict the MRI mode wavevector and mag- 
netic field to be normal to the accretion disk, in 
which case the parallel viscosity has no effect on 
the MRI. This case is the most unstable one in the 
collisional limit, which is the limit in which we are 
mainly interested. 



2. Linear Theory 
2.1. Model 

We consider the two-fluid MHD equations: 
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where 



E = —-v x B + —(-J x B - Vp e ) 



J=^VxB 



ne v c 

P =Pi + Pe- 



The terms representing two-fluid effects are the 
Hall term (J x B/nec) and the electron pressure 
gradient in the definition of E. We assume bary- 
tropic pressure variations of the form 



dp s = TT s dn 



(2) 



for each species s where, for example, T — 5/3 
for an adiabatic equation of state. For l~l we use 
the leading order terms in the Braginskii (1965) 
closure: 

n = it + n gv , 



where 



IT 



|L(|-3»)(S.W.S) (3) 

\b x W • (I + ibb) + (4) 
lu) c i L 

+ [S x W - (I + 366) T | (5) 



where b = B/B, ui ci = eB z0 /miC is the ion cy- 
clotron frequency, v j is the ion collision frequency, 
and the coefficient 770 ~ 0.96 is a factor determined 
by kinetic analysis (Braginskii 1965). The rate-of- 
strain tensor is 

W = X7v + (X7v) T - |lV-u. 

IT is the parallel viscosity term considered by 
Balbus (2004) and Islam & Balbus (2005), which 
represents the lowest-order in d t /vi correction to 
the fluid equations. (It will be shown, however, 
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that the parallel viscosity plays no role in the 
MRI in the configuration we choose to examine.) 
n gv is the gyroviscous force, which represents the 
lowest-order FLR correction to the fluid equations. 
This form of n is appropriate in the limit where 
w c i ^> Vi- Together with the short mean- free-path 
condition (fc||A m / p <C 1, where X m fp ~ Vt/^i) nec- 
essary to justify the fluid closure, this restricts 
the validity of our results to the parameter regime 
where 

k\\v t < v% < Wa. (6) 

In the case of an accretion disk where k\\ > 1/d is 
limited by the disk height d ~ v t /Sl, this validity 
condition becomes f2 <C v% <C io ci . 

2.2. Equilibrium 

We consider a cylindrical equilibrium, using the 
standard coordinates (r, 9, z), and assume axisym- 
metry in the azimuthal coordinate (dg = 0). For 
simplicity, we assume the equilibrium magnetic 
field is of the form B = B za z, and that the equi- 
librium fluid velocity is of the form v a = rCt(r)6. 
Without loss of generality, we orient our coordi- 
nate system so that f2(r) > at the radius of inter- 
est. For such a configuration, radial force balance 
is satisfied when 

9{r)=r^{r) + ^^n'{r)], (7) 

where pi = vt/tOci is the ion Larmor radius, and 
vt = y/Ti/mi is the ion thermal velocity. We ne- 
glect any equilibrium structure in the z-direction. 
In the following derivation, we choose to use equa- 
tion (7) to eliminate g(r) in favor of £l(r) (here- 
after we will drop the explicit dependence of SI on 
r). 

2.3. Local Linear Dispersion Relation 

We consider linear perturbations about this 
equilibrium having scale lengths 1/fc much smaller 
than the equilibrium flow gradient scale length 
L = Q./9, 1 , so that 6 = l/\kL\ < 1. In this limit, 
we may assume that the normal modes of the sys- 
tem are plane waves to lowest order in 5. We re- 
strict our analysis to perturbations cx e 4 ( fcz - wt ) for 
simplicity, as this is the most unstable case in both 
the MHD and collisionless limits (Quataert et al. 
2002). Carrying out the linearization of equa- 
tions (1) yields the following dispersion relation, 



to lowest order in 5: 

= { w 2 + W^^A-C)x 

x {W 4 + a 2 W 2 + a } (8) 

where 

a 2 = -2(2 + R + K) + 

+ 2A{A + R-2A)- H(R + H) 
a Q = [K - A(R + 2H) + 2(R + H)] x 

x[K- A(R + 2H) + H{2 + R)} . 

We define the dimensionless quantities 



W - uo/VL (9a) 

K = (v A k/n) 2 (9b) 

A = Kf3 l (n/uj ci )/A = KF (9c) 

H = K(n/uj ci ) (9d) 

C = {c s k/nf (9e) 

R = rCl'/Cl (9f) 

and characteristic velocities v\ = B 2 /4,imQmi 



and c 2 s = r(T e + T^/rrii. Here fa = 8wp i0 /B 2 
is the ratio of ion thermal pressure to magnetic 
pressure. The dimensionless parameter A mea- 
sures the importance of the gyroviscous force, and 
setting A = is equivalent to omitting V • U gv in 
the ion force equation. Similarly, H measures the 
importance of the Hall term in Ohm's law, and K 
measures the importance of magnetic tension. R 
is the ratio of the radial coordinate to the equi- 
librium flow gradient scale length, and is taken to 
be ~ (9(1). For a Keplerian disk, R = -3/2. For 
convenience we have also defined 

F = ^i — 

4 u C i 

which is the ratio of the gyroviscous force to the 
magnetic tension force. 

Equation (8) contains two uncoupled modes. 
The first factor contains the acoustic mode, which 
may be damped by the parallel viscosity when 
A 7^ 0, and is not of interest here. (While Lo ci jvi is 
generally large, Auj^jvi <~ CSl/vi < Cin the colli- 
sional regime, so the effect on the acoustic mode is 
small.) The second factor, enclosed in braces, con- 
tains the MRI. In the limit where A — > (no FLR 
effects), H — > (no Hall effect), and Cl/i/i — ► 
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Fig. 1. — Contours of the growth rate of the unstable solution to equation (8) are plotted versus Uva/^- 
Left: The growth rate is plotted for various values of F, which measures the importance of FLR effects, in 
the limit where Sl/uj C i <C 1 (no Hall effect). F = is the ideal MHD result. Right: The growth rate is plotted 
for various values of Q/uj c i, which measures the importance of the Hall effect, in the limit where Fcl (no 
FLR effects). Vt/uo ci = is the ideal MHD result. 



(collisional regime), the dispersion relation of Bal- 
bus & Hawley (1998) is recovered. Note that the 
parallel viscosity (oc 7]q) affects only the acous- 
tic mode and not the MRI. Evidently, there is no 
0(Q/ui) correction to the MHD result for the MRI 
when Bq = B z $z, which is in agreement with the 
findings of Sharma et al. (2003) and Islam & Bal- 
bus (2005). Here we are interested in corrections 
to the collisional mode, for which the Bq = B z $z 
case is the most unstable. Extending this analy- 
sis to a more general magnetic field configuration 
substantially complicates the analysis. 

It should also be noted that tu C i is a signed 
quantity since it is proportional to B z , which 
may be positive or negative. Since we have cho- 
sen the coordinate system so that SI is positive, 
sign uj C i = sign B ■ O = sign F. It has been shown 
previously that the effect of the Hall term on the 
MRI depends strongly on the sign relative signs of 
uj ci and SI (Wardle 1999). The effect of the gyro- 
viscous force has a similar dependence. 

The growth rate 7 = Im to of the unstable solu- 
tion to equation (8) is plotted in figure 1 for a Ke- 
plerian rotation profile (R — —3/2). Note that the 
abscissa should be read as a normalized wavenum- 



ber and not a normalized magnetic field strength, 
because F and Sl/w C i are dependent on B. When 
Lad > 0, and hence F > also, both the FLR and 
the Hall effects can be seen to move the most un- 
stable mode to lower wavenumbers, and to reduce 
the value of K at which the MRI is completely 
stabilized. Also, FLR effects increase the growth 
rate of the most unstable mode. When io C i < 0, 
and hence F < 0, both effects are seen instead to 
increase the cutoff value of K all the way to the 
point where modes of any wavelength for which 
this analysis is valid are unstable. Presumably, 
the inclusion of a finite resistivity would damp 
this resonance, as it does in the case of the res- 
onance due to the Hall effect (Balbus & Terquem 
2001). When F or tt/u> c i becomes sufficiently neg- 
ative (F < —2/3 when Q/u) C i — > 0, or SI / cj C i < —2 
when F — > 0), all values of k are suddenly com- 
pletely stabilized. (It has been shown by Balbus 
& Terquem (2001) that this stabilization is less 
sudden when finite resistivity is included.) 

2.4. Instability Criterion 

Applying the Routh-Hurwitz theorem to equa- 
tion (8), we find that the condition for stability of 



4 



an MRI mode is that: 

a 2 < and a > 0. 

This criterion is highly complicated, and for gen- 
eral values of A and H, there may be multiple 
stable and unstable regions in K-sp&ce. 

In the ideal limit, when A — > and H — > 0, the 
instability criterion of Balbus & Hawley (1998), 
K < —2R, is recovered. This limit is well un- 
derstood, and in this case stabilization at high- if 
is due to the effect of magnetic tension. In this 
limit, instability does not exist in flows in which 
the angular velocity increases with radius (R > 0). 

The limit A — > 0, in which the Hall effect is 
dominant over the FLR effects, has also been con- 
sidered before. In this limit, equation (8) reduces 
to the dispersion relation of Balbus & Terquem 
(2001). Since A/H = Pi/4, this limit describes 
accretion disks having /3j -C 1. Formally, the in- 
stability criterion in this case remains somewhat 
complicated because the signs and relative mag- 
nitudes of most of the terms are undetermined 
in general. There is some discussion of the in- 
stability criterion in this case by Wardle (1999) 
and Balbus & Terquem (2001), as well as insight 
into its physical meaning. We will not repeat this 
discussion, except to mention a few interesting 
points. The first is that there may exist some val- 
ues of il/wd for which modes of any wavelength 
are unstable (this is true in the Keplerian case for 
—2 < Q/u> c i < —1/2). Also, some unstable modes 
may be present in disks in which angular velocity 
increases with radius (R > 0), in contrast to the 
ideal result (Balbus & Terquem 2001). 

We are more interested in the opposite limit, 
Pi 3> 1, in which FLR effects are dominant over 
the Hall effect. Taking H — > 0, the dispersion 
relation for the MRI reduces to 

= VF 4 - (11) 
-2 [2 + R + K - A(4 + R- 2A)\ W 2 + 
+ {K- AR + 2R) (K - AR) 

and the criterion for instability is found to be 



For the usual case where R < 0, all modes are 
completely stabilized if F < — l/\R\- As with the 



Hall effect, gyroviscosity allows unstable modes to 
exist when R > 0; in this case, unstable modes 
may exist when F > l/R. In the limit where the 
gyroviscous force dominates the force of magnetic 
tension, F» 1, the instability criterion becomes 
simply 

A < 2 and F > (13) 

or, equivalcntly, (kpi) 2 < 4fl/u> C i, In the case 
where F < 0, there are no unstable modes in this 
limit. 

3. Discussion and Conclusions 

The gyroviscous stress arises from changes to 
the guiding-center drifts due to the gradients of 
magnetohydrodynamic forces (the electric field, in 
this case) across a gyro-orbit. A more complete 
and quantitative discussion of this effect is dis- 
cussed by Kaufman (1960). Due to their larger 
Larmor radius, the ions are more affected by this 
modification that the electrons, leading to the gen- 
eration of currents. Essentially, the stabilization 
(or dcstabilization) of the MRI by FLR effects is 
due to the currents generated by the gyroviscous 
force being out of phase (or in phase) with the 
current of the ideal-MHD MRI eigenmode. 

The effect of the gyroviscous stress becomes rel- 
atively more important to the behavior of the MRI 
as the magnetic field strength is decreased. The 
relative importance of this stress to the magnetic 
tension is A/K = F oc B~ 3 , and the relative im- 
portance of the gyroviscous stress to the Hall effect 
is A/H = f3i/ 4 oc B~ 2 . Note that this means anal- 
yses of the MRI in the "weak-field" limit (in the 
sense that fl /oj C i > 1 ) which did not consider FLR 
effects are valid only for Pi <C 1. If the magnetic 
field is sufficiently weak, the validity condition for 
the Braginskii closure, equation (6), may be vi- 
olated since kvt/oJ c i = \2HF oc in which 
case the above analysis no longer strictly applies. 
Of course, gyroviscosity does not shut off at this 
point; indeed, FLR effects are expected to become 
increasingly powerful as the magnetic field is de- 
creased further, though not in a way that is cor- 
rectly described by equation (4). Therefore it is 
probable that the \ow-K modes which survive at 
the lowest magnetic field strengths covered by this 
analysis will be completely stabilized as the mag- 
netic field decreases further. This is the proper 
resolution to the inconsistency of the ideal MHD 
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result than the MRI remains unstable as B — > 
in the non-dissipative case. Because the MRI is 
not present at arbitrarily low magnetic fields, its 
role in the amplification of primordial astrophysi- 
cal magnetic fields is severely restricted. 

It should also be noted that FLR effects may 
completely stabilize MRI modes having wave- 
lengths much greater than the ion Larmor radius, 
where kpi <C 1. As was shown by Rosenbluth 
et al. (1962), this is possible for "weakly unsta- 
ble" modes like the MRI. Restricting our analy- 
sis to the FLR-dominated limit (F > 1), MRI 
modes are stabilized when (kp^ 2 > 4(D,/u> C i), 
where Q/cJ c j is typically small in astrophysical 
accretion disks. This stabilization may be made 
more important by the fact that in an accretion 
disk the lower bound on k is set by the height of 
the disk d, which may be much smaller than the 
equilibrium flow gradient scale length L. Although 
a proper understanding of this phenomenon must 
take into account the z-stratification of the disk 
equilibrium, we may estimate that the smallest 
wavenumber present in the disk is ~ ir/d. The 
criterion for complete stabilization by FLR effects 
of all MRI modes within an accretion disk at some 
distance from the central mass is then 

^{Pr/df >40/w ci . (14) 

In the typical case where d ~ Vt/Q this inequal- 
ity reduces to u) c i/Q < it 2 / '4. While this is 
not typically satisfied in astrophysical accretion 
disks, it may be satisfied in nascent galaxies with 
weak magnetic fields, or weakly ionized protostel- 
lar disks (Krolik & Zweibel 2006). It is impor- 
tant to recall here that the condition for validity 
of our analysis for wavelengths of this scale re- 
quires SI -c u C i\ whether FLR effects are more or 
less stabilizing than our result would predict out- 
side this range of validity is a matter for further 
research. 

For a concreteness, we consider the example 
of an active galactic nucleus (AGN) in a weak 
magnetic field, for which the conditions for valid- 
ity of the Braginskii equations are met, and in 
which the effect of the gyroviscous force is par- 
ticularly strong. We assume Keplerian rotation 
(R = -3/2), a central mass of 10 7 M©, B = 10 
nG, n = 10 cm -3 , and Tj = 1 eV, at a distance of 
10 16 cm. For this case, the relevant dimensionless 
parameters are F ss 384, H « 3.81 x 10 -4 , and 
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Fig. 2. — The growth rate of the MRI is plot- 
ted versus wavelength for a Keplerian, 1O 7 M0 
active galactic nucleus in a weak magnetic field 
(B = 10 nG), at a distance 10 16 cm. The so- 
lutions to the dispersion relation including (solid 
line), and excluding (broken line) the FLR correc- 
tion are plotted. Gyroviscosity completely stabi- 
lizes the shorter-wavelength modes, and enhances 
the longer-wavelength modes. 

K sa 1.96 x 10~ 5 at k = fl/v t . The growth rates 
of the MRI under these conditions are plotted in 
figure 2. From that figure it can be seen that in 
both cases the FLR correction increases the wave- 
length of the fastest-growing mode from w 10 11 
cm to « 25 x 10 11 cm, and similarly increases the 
cutoff wavelength below which there are no un- 
stable modes. In this example, the weakness of 
the magnetic field is crucial to the importance of 
the FLR effect. If the magnetic field strength is 
raised by two orders of magnitude, the FLR cor- 
rection becomes insignificant. Therefore, "primor- 
dial" accretion disks in weak magnetic fields will 
be most affected by the FLR correction, whereas 
disks in stronger fields (B > 10 pG) may be com- 
pletely unaffected, unless they are unusually hot 
and dense, or have an orbital frequency not much 
smaller than the ion cyclotron frequency. 

The nonlinear saturation of the MRI cannot be 
properly addressed by linear analysis. Because gy- 
roviscosity is not dissipative, and there is no en- 
ergy associated with the gyroviscous term, one 
might expect that the ultimate nonlinear satu- 
rated state should not be affected strongly by the 
gyroviscosity. Ultimately, questions of nonlinear 
saturation should be addressed by numerical sim- 
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illation. 
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